Dynamical Exchanges in Facilitated Models of Supercooled Liquids 
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We investigate statistics of dynamical exchange events in coarse-grained models of supercooled 
liquids in spatial dimensions d = 1, 2, and 3. The models, based upon the concept of dynamical 
facilitation, capture generic features of statistics of exchange times and persistence times. Here, 
distributions for both times are related, and calculated for cases of strong and fragile glass formers 
over a range of temperatures. Exchange time distributions are shown to be particularly sensitive 
to the model parameters and dimensions, and exhibit more structured and richer behavior than 
persistence time distributions. Mean exchange times are shown to be Arrhenius, regardless of models 
and spatial dimensions. Specifically, ~ c~ 2 , with c being the excitation concentration. Different 
dynamical exchange processes are identified and characterized from the underlying trajectories. We 
discuss experimental possibilities to test some of our theoretical findings. 



I. INTRODUCTION 



'lass 



Dynamical arrest of liquids as they approach the ; 
transition is a topic of much current research. 0, 0, li |J| 
This arrest entails such features as non-exponential re- 
laxation and precipitous non- Arrhenius temperature de- 
pendence of transport properties. At its heart lies the 
notion of dynamic heterogeneity. [E IE IE IE IE Eli EH 
Namely, as the temperature decreases towards the glass 
transition, mobility develops spatial inhomogeneity, and 
as time progresses local mobility undergoes dynami- 
cal changes. We have interpreted these phenomena 
0,0,0,0,0] in terms of local excitations of mobility 
in space that pr opag ate in time with facilitated dynam- 
ics. 0, ^E HE l2l| The excitations thus form lines in 
space-time. Dynamic scaling is described naturally in 
terms of the statistics of these lines. 0, 0, |5E 

The work we present in this paper is motivated by re- 
cent experiments that directly detect local fluctuations 
of dynamics in supercooled liquids. 0, IE EH EH 
For example, single molecule rotational experiments 0] 
show that changes in local dynamics are similar to ran- 
dom telegraph noise, where local spatial regions exhibit 
dynamical exchanges between fast and slow dynamics. 
Similar stochastic behavior is also observed in local di- 
electric fluctuation experiments. \J\ In those experiments 
a local spatial region exhibits dynamical exchanges be- 
tween fast and slow dynamics, thus yielding direct con- 
firmation of dynamic heterogeneity. 

Considering these experiments, we establish generic 
features of statistical properties of dynamical exchange 
events and show how these features manifest themselves 
in experimental observables. In the low temperature 
regime, dynamics in supercooled liquids is dominated by 
fluctuations, not the mean. Therefore, it is pertinent to 
study the whole distribution in order to gain insights on 
the microscopic nature of dynamics in supercooled liq- 
uids and glasses. 

We do so in this paper in the context of coarse-grained 
facilitated models. In Sec. HU we define these models. 



In Sec. IIIII we define exchange and persistence times, 
and derive an analytical relationship between the distri- 
butions of exchange and persistence times. In Sec. IIVI 
qualitative features of the exchange time distribution are 
presented. Numerical results of distributions of exchange 
and persistence times for the models are described in 
Sec. In Sec. lVII we present results of numerical simula- 
tions of dynamic bleaching experiments. We conclude in 
Sec. IVIII with discussions of possible experimental mea- 
surements aimed at testing our theoretical predictions. 



II. MODELS OF GLASS FORMERS 

We assume that a kinetically constrained model 0, 
IhH IrsL |5E 0| is obtained through coarse graining over 
a microscopic time scale St (e.g., larger than the molec- 
ular vibrational time scale), and also over a microscopic 
length scale 51 (e.g., larger than the equilibrium corre- 
lation length). The dimensionless Hamiltonian for the 
model is, 



H 



N 



K = o,i). 



(i) 



Here, refers to a state of lattice site i at x», where 
rii = 1 coincides with lattice site being a spatially un- 
jammed region (i. e., carrying mobility), while rij = 
coincides with it being a jammed region (i. e., not car- 
rying mobility). We thus call rii the "mobility field". In 
spatial dimensions d = 1, 2, and 3, the lattice is linear, 
square, and cubic, respectively. The number of sites, N, 
specifies the size of the system. From Eq. (QJ, thermo- 
dynamics is trivial, and the equilibrium concentration of 
defects or excitations is given by 



(rii) 



1 



1 + exp(I/T) 
where T is a reduced temperature. 



(2) 



The dynamics of these models obeys detailed balance 
with local dynamical rules that depend on the configura- 
tion of the lattice site i as well as those of its neighbors. 
Namely, 



t.(+) 



m(t) = o > m{t + St) = 1, 
m(t) = 1 — — > n,(t + St) = 0, 



where 



.(+) _ 



-1/T 



fc, H = MM) 



(3) 
(4) 

(5) 
(6) 



The function fi ({n x }) = /i(ni,ri2, • • • , n^) reflects dy- 
namical facilitation. 

In the Fredrickson- Andersen (FA) model, a change 
of state at site i can occur from t to t + 5t, m(t+St) = 1 — 
rii(t), if at least one of 2d nearest neighbors of the ith site 
is excited. For example, in d = 1 case, 
1 — rii(t) is allowed only if rij(t) = 1 where xj = x, ± ^i. 
In contrast, in the c£ = 1 East model case, [2(| one needs 
an excited neighbor in specific directions in order to have 
a state change at site i. For instance, rij(t) = 1 such that 

X j — X i 5 1/ • 

In d = 3 cases, rij(t) at Xj = (x, y, z) can make a change 
if at least there is one nearest neighbor site j such that 
rij(t) — 1 at Xj = (xiSl, y, z) or (x, y±5l, z) or (x, y, z± 
<5^) (FA model) or Xj = (x — Sl,y,z) or (x, y — <5Z,z) or 
(x, y, z — SI) (East model). The facilitation function for 
the FA and East model, fi ({n x }), is therefore given by, 



[/i({nx»] 



East 



i-IK 1 



(7) 
(8) 



where u; is a unit vector in the Ith dimension. In the clas- 
sification scheme of kinetically constrained models given 
in Ref. l25l these kinetic constraints correspond to one- 
spin facilitation rule in d dimension. 




FIG. 1: Exchange and persistence times are shown in the 
trajectory of d = 1 FA model at T — 0.8. Shaded regions 
represents parts of space-time with excitations, while white 
regions represent parts with no excitations. 



that the decoupling of the translational diffusion from 
the relaxation times in supercooled liquids near the glass 
transition can be explained from the distribution of ex- 
change times. 0, Persistence and exchange times 
are illustrated in the trajectory of d = 1 FA model in 
Fig.ID 

Exchange and persistence times are different statisti- 
cal measures of the same trajectories. There is a general 
relation between exchange and persistence time distribu- 
tions, independent of models. To derive the relation, we 
define the variable 



Pi(t;l,r) = 



t-5t 



n ^■( r + < ') 



[i-m(T + t)}. (9) 



It is the density of excitations at time r (i.e., m = 1) 
that persist until time t + r. A similar expression for the 
density of unexcited sites that persist for this time frame, 
Pi(t; 0, r), is given by Eq. © with changed to 1 — n,-, 
and 1 — Hi changed to n^. 
Similarly, 



Xi(t;l,T)= I[l-^(r-<5i)] 
x [l-mir + t)], 



t-5t 



n n *( T+t ') 



t'=0 



(10) 



III. EXCHANGE AND PERSISTENCE TIMES 

Relaxation time scales of the glassy systems can be 
described by a distribution of persistence times, t p , the 
time at which a local region changes its state for the first 
time once the trajectory has started at time zero. In the 
d = 1 FA model, for instance, the mean persistence time 
exhibits Arrhenius behavior at low temperatures. |T^ 

A statistical measure which is useful in studying 
changes in dynamics of local environments is the distri- 
bution of exchange times, t x , the time duration of given 
states of mobility fields. For example, it has been shown 



is the density of excitations that exchange over the time 
interval between r and t + r. Clearly, 



Xi(t;l,T) 



1 

Si 



[Pi(t;l,T)-Pi(t + 6t;l,T-5t)]. (11) 



Averaging and taking St — > + therefore yields 

dV(t; 1) 



X(t;l) = - 



dt 



(12) 



where X(t; 1) = (A^;1,t)) and V{t;l) = (Pi(t;l,r)). 
The normalized probability densities of exchange and 
persistence times for an excitation, p x (t] 1) and p p (t; 1), 



respectively, are proportional to X(t; 1) and V(t; 1), re- 
spectively. Accounting for normalization constants there- 
fore gives 



Pp 



(t;l) = £W) 



(*x(l)> 



(13) 



where (t x (l)) = J °° di tp x (t; 1) is the mean exchange time 
for an excitation. 

One can also proceed the same route for persistence 
and exchange of no excitation by replacing m with 1 — rii 
in Eqs. © and (QDJ to find 



Pp(*; o) 



(*x(0)> ' 



(14) 



The overall exchange and persistence time probability 
densities, p x (t), and p p (t), are given by 



Px(t) = 2 M*;0)+fe(*;i)]» 



(15) 



Pp (t) = (l-c)p p {t;0)+c Pp (t;l). (16) 

Using Eqs. i|13fl and (|14|l with the condition of detailed 
balance, 



we find 



P P (t) 



(*x(l)) _ C 

<*x(0)> 1-c' 



(17) 



(18) 



where (t x ) = |((£ x (0)) + (t x (l))). For a Poisson process, 
which is valid for unconstrained dynamics, p x (t; n) = 
r(n)~ 1 exp(— t/r{n)). In that case, the exchange and 
persistence time distributions become identical to each 
other. 

Moments of the two distributions are related to each 
other from Eq. 118J1 through integration by parts. In 
particular, 



((t P y 



(m + l)!(t x )' 



(19) 



The moments of the persistence time distributions are 
always greater than those of the exchange time distribu- 
tions when the distributions are broader than Poissonian. 
Correlations between exchange events are described by 
correlation functions of the type 

C l3 {t,t'\n,T,n',T') = (8Xi(t;n,T)5Xj(l!;n'S)), (20) 

where SXi(t; n, r) = Xi(t; n, r) — X(t; n). 



IV. ANALYSIS OF EXCHANGE TIME 
DISTRIBUTIONS 

We have calculated exchange and persistence time dis- 
tributions for the FA and East models by performing 




2 2 4 6 
lo g,nt 



FIG. 2: Exchange and persistence time distributions for d = 1 
FA and East models. Here, P x (log 10 t) = tp x (i) ln(10), and 
P p (log 10 t) is similarly related to p P (t). Exchange time dis- 
tribution from simulation (filled circles) , persistence time dis- 
tribution from simulation (open circles), and persistence time 
distribution predicted from Eq. 11811 (solid lines), (a) d = 1 
FA model at T = 1, (b) d = 1 FA model at T = 0.5, (c) d = 1 
East model at T — 1, and (d) d = 1 East model at T = 0.5. 



Monte Carlo simulations. For the purpose of numeri- 
cal efficiency, we have used the continuous time Monte 
Carlo algorithm. [27il28| To sample long exchange times 
(i.e., log 10 t x > 1), we have used N = 100c _1 . To sam- 
ple short exchange times (i.e., log 10 t x < 1), we have 
used N = 10 5 c _1 . The two distributions obtained were 
matched at log 10 t x — 1. Simulations were performed for 
total times T — 500 tq, with ra being the relaxation 
time of the model. Averages were performed on about 
100 independent trajectories in each case. 

As an illustration, we show in Fig. [2] exchange and 
persistence time distributions obtained from numerical 
simulations, and compare those results with the predic- 
tion of Eq. l(TH|l . The comparison verifies Eq. ljTB|) . The 
peak and principal statistical weight in the persistence 
time distribution occurs at a longer time than that for 
the exchange time distribution. 



A. Moments of Exchange and Persistence Times 

In many experiments, measurements of moments of 
fluctuating quantities are more easily made than mea- 
surements of the whole distributions. Studies of mo- 
ments will yield information on underlying dynamics. 
p9T l30| We have studied temperature dependence of the 
moments of exchange and persistence time distributions, 
((t x ) m ) and ((t p ) m ) (m=l, 2, and 3), respectively, of the 
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FIG. 3: The first (circles), second (squares), and third (dia- 
monds) moments of exchange (open symbols) and persistence 
(shaded symbols) time distributions of d = 1,2, and 3 FA 
models (left panels) and East models (right panels) as func- 
tions of excitation concentration, c. 



FA and East models. Figure 01 shows these moments for 
dimensionality d—l, 2, and 3. 

Moments of the exchange and persistence times re- 
veal different temperature dependence for the FA and 
the East models. In the FA models, all the moments 
of exchange and persistence times increase in an Arrhe- 
nius fashion, (t m ) ~ c _Qm , as the temperature (or the 
excitation concentration) decreases. In the East mod- 
els only the mean exchange times are Arrhenius. Higher 
moments of exchange time distribution and all the mo- 
ments of persistence time distribution in the East models 
are super- Arrhenius. Further, the mean persistence time 
is larger than the mean exchange time in all cases. Table 
H] lists scaling exponents of the moments for the exchange 
and persistence times for Arrhenius cases. 

We find that the mean exchange time increases approx- 
imately as 



(21) 



in the low temperature limit in all cases, regardless of the 
models and dimensionality. To understand this result, we 
use Eqs. Q and © to establish that the average value 



FA models 




m = 1 


m = 2 


m — 3 


d= 1 


2.0(3.2) 


5.2(6.4) 


8.4(9.7) 


d= 2 


1.9(2.3) 


4.2(4.7) 


6.6(7.0) 


d = 3 


1.9(2.1) 


4.0(4.2) 


6.1(6.3) 



East models 




m — 1 


d= 1 


2.0 


d= 2 


2.0 


d = 3 


1.9 



TABLE I: Scaling exponents for fitting various moments of ex- 
change and persistence time distributions to Arrhenius form, 
(t m ) ~ c~ a ™ are shown for d = 1, 2, and 3 FA and East 
models. Numbers in the parenthesis are those for persistence 
time distributions. Fitting of the moments to the Arrhenius 
form was done for data corresponding to T < 1. 



of facilitation function, (fi), is 

(/,) = !-(!- c) ad , 



(22) 



where a = 2 and 1 for FA and East models, respectively. 
As such, the mean exchange times for n = and n = 1 
states obtained from flip rates in Eqs. 10 and © are 



<*x(0)> 

<*x(l)> 



l/(A ; | +) )=cxp(l/T)/(/ i ), 



(-)\ 



!/</<>, 



(23) 
(24) 



and the overall mean exchange time, (i x ) = ^((i x (0)) 
(**(!))), ^ 



1 



1 



1 



as c 



0, 



2c(/i) 2c[l - (1 - c) ad ] 2acfc 2 ' 

(25) 

which shows the observed scaling of (i x ) ~ c 2 in the 
low temperature limit with a proportionality constant 
that depends on the dimensionality and model. Figure 01 
shows that both the scaling and the proportionality con- 
stant predicted within Eq. i12">l) agree with the simulation 
results. 

This behavior of the mean exchange time versus the 
excitation concentration has at least two experimental 
implications. First, experimental measurements that are 
insensitive to fluctuations of exchange times will not dis- 
tinguish between fragile and strong glass formers. In or- 
der to distinguish, measurements of higher moments, or 
perhaps, the whole distributions will be required. The 
distribution could possibly be observed through single 
molecule experiments. [T3, 122, y(| Second, the universal 
temperature dependence of the mean exchange time can 
form a basis of experimental measurements of excitation 
concentrations in real systems. The specific c -2 scaling 
holds to the extent that experimental systems can be 
modeled with single spin-flip facilitation rules. 

In contrast to the behavior of the mean exchange 
time, the temperature dependence of mean persistence 
time, (t p ), depends on the model and the dimensional- 
ity. Relaxation times of the models are given by their 
mean persistence times, or from Eq. I|19|) . the second mo- 
ment of the exchange times. Temperature dependence of 
moments of exchange and persistence times shows that 
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1000 



FIG. 4: Comparison of mean-field approximation, Eq. 125H 
(dashed line), and simulation results of mean exchange times 
in d = 1, 2, and 3 FA and East models, a — 2 for FA models 
and a — 1 for East models. 



relaxation times are Arrhenius for the FA models and 
super- Arrhenius for the East models. 

As the dimensionality of the model increases, differ- 
ences between the moments of exchange and persistence 
times become smaller in the FA model. The model be- 
comes increasingly mean-field like in higher dimensions. 
No such convergence is found in the East model for all 
dimensions investigated. The East models have weak di- 
mensional dependence due to the quasi-one dimensional 
nature of directional persistence. 



B. Stokes—Einstein Breakdown 

For the FA models, which correspond to strong glass 
formers, the translational diffusion constant scales as 0] 



D 



FA 



(**) 



(26) 



This result is true in all dimensions, ft follows from 
Eq. (|21[) and the statement that adjacent random walk 
steps coincide with adjacent exchange events. For the 
East models, however, the result does not hold because 
adjacent exchange events are correlated. In particular, 
for the hierarchical dynamics of the East models, large 
bubbles in space-time are covered with smaller bubbles, 
which in turn are covered by yet smaller bubbles, and so 
on. In that case, higher moments of exchange time 
distributions exhibit different temperature dependence 
from that of the mean. This behavior is illustrated in 



Fig. 13 for the East model, which corresponds to a fragile 
liquid. Correlations between successive exchange events 
can be investigated by using correlation function given 
in Eq. I|20[l . The effects are responsible for non-trivial 
scaling relation between the translational diffusion and 
persistence times in fragile liquids. 0, |35j |3(| lU ^ 
Indeed, numerical results for the fractional Stokes- 
Einstein relationship, 



(27) 



give £ « 2/3, 2/2.3, 2/2.1 for dimensions d — 1, 2, 3 in the 
case of the FA models, while £ rj 0.7 — 0.8 (very weakly 
dependent on d) in the case of the East models. l31| 



C. Classifications of Exchange Events 

We show in Fig.[S]representative exchange time distri- 
butions. At short exchange times, log 10 t x < — 1, the dis- 
tribution P x (log 10 frx) follows a power-law behavior with 
a slope 1 in all cases, P x (log 10 t x ) ~ log 10 t x . This is be- 
cause P x (log 10 t x ) ~ t x p K (t x ), and Poissonian statistics is 
obeyed, p x (t x ) ~ exp(— t x ), for short times. 

Multiple peak structures develop in the exchange time 
distributions at low temperature. In order to investigate 
the structure of exchange time distributions in detail, we 
define mobility index, m,i(t), that depends upon whether 
the lattice site is filled with an excitation or not and 
also whether it is in an mobile or immobile configuration. 
Specifically, 



mi(t)=4-2m(t)-fi(Jt), 



(28) 



where fi(t) = 1 when the lattice site i is in a mobile 
configuration at time t, and it is zero otherwise. As such, 
rrii{t) is 1 when site i is excited and mobile, 2 when it is 
excited and immobile, 3 when it unexcited and mobile, 
and 4 when it is unexcited and immobile. 

During an exchange time of the lattice site i, rij does 
not change, while its nearest neighbors may make flipping 
events, thus changing rtii over time. In order to monitor 
changes in the local mobility during an exchange time, 
it is useful to introduce an averaged mobility index for a 
single exchange event, 

mi{t s ) = i / dt' mi (t s + t'), (29) 
t x Jo 

where t s is the start time of an exchange event for site 
i, and the zero time in this formula is the time at which 
the exchange event begins. See Fig. ^ If there has not 
been any change in the local configuration during an ex- 
change time, the lattice site remains always mobile, and 
the averaged mobility index will be either m, = 1 (for 
rii = 1) or m, = 3 (for = 0). However, when the local 
mobility changes during the exchange time period due to 
changes of the states in nearest neighbors, the averaged 
mobility index will be either 1 < < 2 (for rij = 1) or 
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FIG. 5: Decompositions of exchange time distributions of FA 
models (left panels) and East models (right panels) for d — 
1,2, and 3 at T — 0.3. Four sub-distributions - case 1 (circles), 
case 2 (squares), case 3 (diamonds), and case 4 (triangles) - 
add up to the full distribution (solid line). Positions of three 
peaks in P(log 10 i x ) are shown for d — 3 FA model cases. 



3 < rrii < 4 (for rij = 0). We define four different cases 
of the averaged mobility. Namely, 



case 1 
case 2 
case 3 
case 4 



nil = 1 {rii = 1 and /j = 1) 

I < rrii < 2 (rii = 1 an( i /i changes) 

m, = 3 (rii = and fi = 1) 

3 < m t < 4 (rii = and /j changes). 



(30) 



Examples of exchange times that belong to each case are 
shown in the trajectory picture of d = 1 FA model in 
Fig. 

To illustrate the correlation between exchange times 
and environments, Fig. [3] shows exchange time distribu- 
tions and their sub-distributions, each corresponding to 
one of the cases in Eq. (|30|l , for d = 1, 2, and 3 FA and 
East models at T = 0.3. There exist multiple peaks in 
the distributions of exchange times at low temperature 
in all the cases shown in Fig. [SJ Positions of peaks in 
the distribution, t*, and t%, are shown in the case 
of d = 3 FA model. Temperature dependence of these 
peaks for d = 1, 2, and 3 FA and East models are pre- 




FIG. 6: Examples of exchange events that belong to each case 
of Ea. 13UH are illustrated in the traiectorv of d — 1 FA model. 
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FIG. 7: Positions of peaks in the exchange time distributions 
are shown for various excitation concentrations. 



sented in Fig. and in all the cases, each peak position 
is Arrhcnius, 



t: 



(31) 



where n = Q, 1, and 2. 

To reveal the physical origins of each peak in the ex- 
change time distributions, the temperature dependence 
of each peak in the exchange time distribution is shown 



in Fig. Cases 1 and 2 will involve disappearances of ex- 
citations, while cases 3 and 4 appearances of excitations 
in the trajectory space. The first peak in the exchange 
time distribution corresponds to fast processes of mobile 
n = 1 states embedded in the excitation line. In this case, 
the excited state will quickly de-excite in time ig ~ c °j 
and the position of the first peak in the exchange time dis- 
tribution does not depend on the temperature as shown 
in Fig. [7] This process will have the averaged mobility 
index, m = 1 , corresponding to case 1 in Eq. i|3U|) . 

The mobile region located at the boundary of excita- 
tion line will have a probability of becoming immobile 
due to changes in its nearest neighbors and becoming 
mobile later. This process corresponds to case 2 with 
the averaged mobility index 1 < m < 2. It is respon- 
sible for fluctuations of thicknesses of excitation lines in 
trajectory space, and it gives rise to the second peak in 
the exchange time distribution, with temperature depen- 
dence t* ~ c _1 . 

Dynamical exchange processes corresponding to case 3 
also contribute to the second peak in the distribution. In 
this case, n = state next to an excitation line becomes 
excited in time t* ~ c . This process involves a creation 
of an excitation next to a pre-existing excitation line. In 
terms of the geometry of the trajectory space, exchange 
events corresponding to cases 2 and 3 are responsible for 
bending of excitation lines. 

Finally, the exchange events that have the longest 
times, resulting in the third peak, will correspond to the 
up-flip event of the n = states far away from the exci- 
tation lines. This case has the averaged mobility index 
of case 4. In d = 1 FA model, for example, it corre- 
sponds to a region deep inside the bubble structure in 
the trajectory space. 



Exchange times Persistence times 




FIG. 8: Distributions of the logarithm of exchange and per- 
sistence times of FA model for d — 1,2, and 3 at various 
temperatures. From left to the right: T = 10, 1, 0.55, 0.4, 
0.3, 0.22, 0.18 for d = 1 case, and T = 10, 1, 0.55, 0.4, 0.3, 
0.22, 0.18, 0.16 for d = 2 and d = 3 cases. 



V. DISTRIBUTIONS OF EXCHANGE AND 
PERSISTENCE TIMES 

A. Strong Glass Former Models 

Figure[S]compares the exchange time distributions (left 
panels) and persistence time distributions (right panels) 
for d = 1,2, and 3 FA models at various temperatures. 
Exchange time distributions shown in the left panel de- 
pend strongly on the temperature. In the high tempera- 
ture regime, the exchange time distribution has a single 
peak because there is no clear distinction between dif- 
ferent exchange events due to the mean-field nature of 
the dynamics in the high temperature regime. However, 
as temperature decreases, the mean-field picture is no 
longer valid, excitation lines dominate trajectory space, 
and different kinds of dynamical exchange events mani- 
fest themselves in the distribution as distinct peaks. This 
fact is evident from the multiple peak structures in the 
exchange time distribution in all dimensions in the FA 
model. In d = 3 FA model, for instance, distinct triplet 
structure is observed when T = 0.4 and below. Each peak 



corresponding to different dynamical exchange event ex- 
hibits different behaviors of temperature dependence as 
presented in Fig. [7| 

As the dimensionality increases, the exchange time dis- 
tribution shows more and more well-separated triplet 
peak structure in the low temperature regime. The rela- 
tive weight of the third peak increases as the dimension- 
ality increases. 

The persistence time distributions shown in the right 
panel of Fig.[S]are less structured than the exchange time 
distributions, as to be expected from Eq. ljl"%|) . As the 
temperature decreases, the single peak in the persistence 
time distribution moves towards the long time region. 
Physically, persistence time distributions are statistical 
measures of duration of a state at an arbitrary space-time 
point in the trajectory space. Thus, relative to exchange 
time distributions, persistence time distributions empha- 
size processes occurring over longest times. The greater 
structure in the exchange time distributions, relative to 
the persistence time distributions, is due to the former 
giving equal weights to events, irrespective of their dura- 
tions. 
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FIG. 9: Distributions of the logarithm of exchange and per- 
sistence times of East model for d — 1, 2, and 3 at various 
temperatures. From left to the right: T = 10, 1, 0.6, 0.5, 0.4, 
and 0.3 for d = 1 case, T=10, 1, 0.5, 0.4, 0.3, 0.25, and 0.22 
for d = 2 case, and T=10, 1, 0.5, 0.4, 0.3, 0.25, 0.22, and 0.2 
for d = 3 case. 



B. Fragile Glass Former Models 

Exchange and persistence time distributions for the 
East model cases are shown in Fig. |5J In the high tem- 
perature regime, the exchange time distributions have a 
single peak as in strong liquids or FA models. However, 
in the low temperature regime, in contrast with its be- 
havior in the FA model, they show two dominant doublet 
peaks in all dimensions. We also notice that, albeit small, 
broad tails develop in the wings of the distributions at 
low temperatures. This feature is especially evident in 
Fig. ^3 where the exchange time distributions of the FA 
and East models are compared in log-log scale. In the 
case of persistence time distributions shown in the right 
panel of Fig. |5J single peaks develop into broad distribu- 
tions with shoulders as temperature decreases. 

We compare exchange time distributions in FA and 
East models in Fig.^3 In d — 1 case, exchange time dis- 
tributions of both FA and East models show long time 
tails in the low temperature regime in common. Addi- 
tionally, the distributions in the East model exhibit log- 
arithmic oscillations. In higher dimensions, d = 2 and 
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FIG. 10: Exchange time distributions of FA models (left 
panels) and East models (right panels) in log scales. Tem- 
peratures are the same as those given in the captions to 
Figs. H and El A power-law corresponding to p x (ix) ~ t x a 
or P x (log 10 f x ) ~ t x p x (i x ) ~ i]r a is shown as a dashed line 
with a — 0.52, and a stretched exponential corresponding to 
Px(tx) ~ e~'* x ' where j3 = 0.38 is shown as a dot-dashed line 
for d — 1 FA model. In d = 2 and 3 FA models, exponential 
fits, p x (tx) ~ e~ tx ^ T , are given as dashed lines for the lowest 
temperature of each case. 



3, the exchange time distributions in FA models exhibit 
distinct multiplet peak structures due to their mean-field 
natures in higher dimensions. However, the distributions 
in the East models remain similar to those in d = 1 cases, 
exhibiting long time tails and superimposed logarithmic 
oscillations. These are due to the hierarchical dynamics 
that arises from directed facilitation in the East models 
at all dimensions. 

In the d = 1 FA case, due to the formations of bubble- 
like structures in the trajectory space, the exchange time 
distributions decay with slowly decaying long time tails. 
It turns out that in d = 1 FA model case the exchange 
time distribution first decays as a power-law (dashed line) 
then turns into a stretched exponential (dot -dashed line) 
as shown in Fig. However, in d = 2 and d = 3 FA 
models, they decay as an exponential at the long time re- 
gion shown as dashed lines. This indicates that in higher 



FIG. 11: Persistence time PDF's at a high (T = 10) and a 
low (T = 0.3) temperature at equilibrium (solid lines) and 
right after the bleaching (dotted lines). Each dashed line cor- 
responds to the distribution created by bleaching efficiency 
function, 4>{i) with different values of r, given in the inset 
of Fig. 1111 . Inset: Bleaching efficiency function used in the 
simulation, Eq. l-'i'it . with a = 0.08. From left to the right, r 
progressively increases from r = 10~ 6 to t = 10 6 by a factor 
of 100. 



dimensions the FA models are less "glassy" than d = 1 
case. In higher dimensions, chances of forming bubble- 
like structures will be much smaller than in d = 1 case. 
In the East model cases, similar behaviors are observed 
in all dimensions. Namely, the exchange time distribu- 
tions decay with multiple transitions superimposed with 
stretched exponentials in the low temperature regime. 



VI. SIMULATIONS OF BLEACHING 
EXPERIMENTS 

Experiments have been performed in order to mea- 
sure the timescale of the approach of non-equilibrium 
distributions of local relaxation times to equilibrium. 
p^. l3"3 . l33T,l34j | The timescale at which the non-equilibrium 
distribution approaches to equilibrium will be called the 
recovery time. Based on the d = 1 East model, we present 
results of numerical simulations of one of such experi- 
ments - dynamical hole burning experiments or bleaching 
experiments. 

In numerical simulations of the bleaching experiment, 
the bleaching efficiency function, <f>(t), is introduced as a 
probability of bleaching out the region that has a local 
persistence time of t. We choose the following form as 



cf>(t) = min [l,(r/t)°] 



(32) 



where r and a can be varied. We choose r = 1, and a — 
0.08 motivated by experimental results in Ref.|32- After a 
lattice site is chosen to be "bleached out" in simulations, 



FIG. 12: Mean persistence times right after the bleaching ex- 
periment, tobs(0), are shown for various bleaching depths and 
temperatures. £ o b s (0) increases more abruptly as a function 
of the bleach depth at a lower temperature. 



persistence times of that site are not collected at later 
times as time goes on. 

Numerical values of the bleaching efficiency functions 
are shown in the inset of Fig. ^2 The bleaching depth, B, 
is defined as the fraction of lattice sites that are bleached 
out by applying the bleaching efficiency function, and is 
given by 



B = 



dtcb(t)p p (t). 



(33) 



Probability distributions of persistence times after the 
bleaching are shown in Fig. ^2 for two different values of 
the temperature and various values of r. 

Immediately after the bleaching, the system is out of 
equilibrium, with a non-equilibrium persistence time dis- 
tribution given (to within a normalization constant) by 



Pne(t) = 



p p (t)(l-<f>(t)) 



1 - B 



(34) 



After waiting a time i w , this distribution relaxes to 
p(i;i w ), where p(t;0) is p ne (t) and p(t; oo) is the equi- 
librium p p (t). tobs(iw) is the first moment of p(t; t w ), 



tohs (^w) 



dt t pit; t v 



(35) 



and ^obs(O) > (t p ). iobs(0) at various bleaching depths 
and temperatures are shown in Fig. 1121 

The relaxation dynamics is monitored by following a 
time evolution of the non-equilibrium persistence time 
distribution obtained from un-bleached lattice sites for 
various waiting times, t w . Figure H31 shows the relax- 
ation dynamics of the initial non-equilibrium distribu- 
tion of persistence times for d = 1 East model at a low 
temperature. 




FIG. 13: Relaxation dynamics of the distributions of persis- 
tence times to equilibrium after the bleaching of d = 1 East 
model at T = 0.5. r = 1 and a = 0.08 were chosen in <p(t). 



The timescale of the relaxation process of the non- 
equilibrium distribution can be estimated by several 
means. For example, one can look at the waiting time 
dependence of moments of the distributions. Also, it is 
possible to measure how far a non-equilibrium distribu- 
tion is separated from equilibrium by defining a distance 
between two distributions by 

/>oo 

A(i w )= / dt\p(t;t w )-p p (t)\. (36) 
Jo 

A(£ w ) and the first three moments of the non-equilibrium 
distribution are shown in Fig. 1141 for different values of 
the waiting time. From these calculations, the timescale 
of a recovery to equilibrium is estimated to be on the or- 
der of the equilibrium mean persistence time, (t p ). This 
predication is at variance with results from bleaching ex- 
periments, 0, l3^| but is is consistent with results from 
NMR experiments. (3^, 0, The reason(s) for this 
discrepancy remains to be clarified. 



VII. DISCUSSION 

In this paper, we have provided detailed illustrations 
of the generally large differences between distributions of 
persistence times and exchange times. The distributions 
are the same only when these distributions are Poisso- 
nian. Thus, the differences are manifestations of dynamic 
correlations, correlations that exist due to dynamic het- 
erogeneity. Despite the differences between the two dis- 
tributions, the behavior of one determines the other, as 
Eq. iJTSJ shows. 

The first moment of the exchange time distribution 
is Arrhenius, independent of model and dimensional- 
ity. This result can be important if an experimental 




FIG. 14: Recovery of the non-equilibrium distribution, 
p(i;t w ) to equilibrium is shown by the decay of A(£ w ) (up- 
per panel) and-Rm(iw) = ((iobs(£w)) m )/((£ P ) m ) (lower panel) . 
Timescale of the recovery to equilibrium is on the order of the 
equilibrium mean persistence time itself. The model parame- 
ters are the same as those used in Fig. ED 



means is available to measure this moment. In partic- 
ular, Eq. 21 shows that the measurement of the mean 
exchange time determines the average mobility concen- 
tration, c, the central control parameter in our descrip- 
tion of glassy dynamics. [13|, EH, [H, EH E3 This can be 
accomplished with single-molecule experiments in which 
exchange times for a probe molecule passing between fast 
and slow environments are measured directly. [12] The 
determination of c in principle enables parameter-free 
tests of scaling relations we have predicted for transport 
properties. 

Significant experimental manifestations of the dif- 
ferences between exchange-time distributions and 
persistence-time distributions are decoupling of different 
measures of dynamics. The breakdown in Stokes- 
Einstein relations discussed in this paper is one example. 
Another is the existence of the Fickian crossover, 
marking the coarse-graining length-scale beyond which 
random- walk diffusion is correct. [TtI ] At large wave- 
vectors, dynamics is dominated by persistence and 
molecular motion is not Fickian. At small wave-vectors, 
dynamics is dominated by exchange, and molecular 
motion is Fickian. 

Thus, there are many different relaxation times for 
slow dynamics in a glass former. Those connected to 
the mean exchange time will be Arrhenius, even for frag- 
ile glass formers. Non- Arrhenius behavior (i.e., fragility), 
when it exists, is described by the second and higher mo- 
ments of the exchange-time distribution. This feature of 
dynamics in glass formers is general, not a consequence 
of a particular model. It follows from the presence of 
dynamic heterogeneity and its consequential fluctuation 
dominance of dynamics. 
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